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The description of the properties of liquid Helium is a challenge for any microscopic 
many— body theory. In this context, we study the ground state and the excitation 
spectrum of one ^He impurity in liquid '^He at T = with the aim of illustrating 
the power of the correlated basis function formalism in describing heavily corre- 
lated systems. The strong interatomic interaction and the large density require 
the theory to be pushed to a high degree of sophistication. A many— body corre- 
lation operator containing explicit two- and thre-particle correlation functions is 
needed to obtain a realistic ground state wave function, whereas a perturbative 
expansion including up to two phonon correlated states must be enforced to study 
the impurity excitation energies. The theory describes accurately the experimen- 
tal spectrum along all the available momentum range. As empirically shown by 
the experiments, a marked deviation from the quadratic Landau-Pomeranchuck 
behavior is found and the momentum dependent effective mass of the impurity 
increases of ~ 50 % at ij ~ 1.7 with respect to its q = value. Although 

the main emphasis is given to the Correlated Basis Function theory, we present 
also comparisons with other methods, as diffusion Monte Carlo, variational Monte 
Carlo with shadow wave functions and time dependent correlations. 

1 Introduction 

Atomic Helium fluids are an endless source of physical motivations for both theo- 
rists and experimentalists. We have so far accumulated a huge amount of experi- 
mental information and the activity is continuously pushed further to, explore new 
situations. An cpcample is the physical realization of almost two- tl (films) and 
one-dimensional (nanotubes) systems and the consequent possibility to study the 
dependence on dimensionality of correlation effects. Helium clusters are another 
extremely promising field and the present experimental capabilities make pcesible 
to investigate the minimum number of atoms needed to have superfluidity u. For 
theoreticians, Helium liquids can be considered as excellent laboratories to test 
many-body theories. In fact, the interaction is simple and depends only on the dis- 
tance between the atoms. In addition, the effects of having different kinds of quan- 
tum statistics may be studied, since here the quantum behaviour has macroscopic 
manifestations. Actually, the fact that Helium remains liquid at zero temperature 
is a consequence of the large zero point motion of the atoms in the fluid and can 
be considered as a macroscopic quantum effect. p. 

In spite of the enormous progress in the last two decades □ we have not yet 
achieved a completely satisfying explanation of some experimental facts. Let's just 
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mention that a full microscopic description of superfluidity in liquid '^He is still 
missing. It is not possible, in a single presentation, to summarize the present 
status of the whole Helium physics field, so we will exploit the fact that in the 
last two years several microscopic many-body theories have been used to study the 
ground and excited states of one '^He impurity in liquid "^He at zero temperature 
and use this system to illustrate the power and peculiarities oLthese techniques. 
The main emphasis will be given to Correlated Basis Fupction U'Q but comparisons 
with the results obtained with Shadow Wave jFunctionEl (SWF), Time Dependent 
Correlations ^(TDC) or Difusion Monte Carloy(DMC) wiU be given. 

The theoretical study of one ^He impurity in atomic liquid ^He is very helpful in 
understanding the properties of '^He-'*He mixtures. "^He and ''He are isotopes which 
follow Fermi and Bose statistics, respectively. They can form a liquid mixture 
which remains stable at zero temperature, with a maximum solubility of the ^He 
component at zero pressure of ^^6.6%. Both types of statistics coexist in the mixture 
with the consequence that the excitation spectrum is particularly rich. There are 
two types of excitations, whose lowest energy one corresponds to a particle-hole 
band associated to ^He quasi-particles. The quasi-particles have a single particle 
spectrum characterized by an effective mass, mg, mainly due to the interaction with 
the '^He atoms. A little higher in energy lies the collective phonon-roton branch, 
associated to ''He and little affected by the presence of the '^He component. There 
is a marginal low momentum admixture between the two branches and they have 
been clearly separated in recent inelastic neutron scattering experiments at low 
momentum transferlllj. 

As the '^He concentration, X3, in the mixture is small, it is justified to study the 
X3 — »0 limit, or the impurity problem. The analysis of the impurity behavior not 
only provides clues to understand the mixture itself but also gives useful information 
on pure "^He liquid. As we will see, the impurity can be used as a theoretical probe 
to examine the kinetic and potential energies of the bulk Ej system. 

The bulk properties of helium liquid are well measured 113, ^ajid a complete 
thermodynamic information is also available for ^He-^He mixtures i3. In particular, 
the chemical potential of the ^He impurity at ^He saturation density (poC^He) = 
0.02186A^'^) and zero temperature is /i3 = —2.78 K, to be compared with the 
chemical potential of pure ''He , (14 = —7.17 K, and with the one of liquid ^He, 
/i3 = —2.5 K at its own saturation density, po{^He) = 0.0|1£3A^^. 

The recent inelastic neutron scattering experiments £3 at low momenta and 
low concentrations have given access to the impurity excitation spectrum. The 
experiments shcpa a significant deviation from the quadratic Landau-Pomeranchuck 
(LP) spectrum OJ, eLp{q) = Ti'^q^/2m\. The measured spectrum is well described 
in a modified LP (MLP) form 



27713 1+79^ 



eMLp{q) = 7^ 1 , 9 ■ (1) 



Although there are still some uncertainties in the analysis of the data 0, the esti- 
mated values of the MLP parameters, at_E = 1.6 bar and X3 ~ 0.05 are ml y-2-'2 ma 



and 7 ^ 0.13 A^. A recent analysis Ej of specific heat measurements E3 gives 
a slightly different zero momentum effective mass of 7713 = 2.18 m^, taking into 
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account the corrections due to the '^He finite concentration in the mixture. 

The purpose of this contribution is to report the recent progress in the micro- 
scopic description of the impurity system. This wiU give us the opportunity to 
comment also on the present status of the description of the "^He hquid and enhght 
some aspects of the dynamics of ^He-'*He mixtures. 



2 Static properties 



For a quantum microscopic description we start from an empirical Hamiltonian, 
defined in terms of the masses of the atoms and of their mutual interactions. In 
the case of one impurity, the Hamiltonian reads 

^ - - E + E ^(1 ^ rf I) - + E ^(1 r(3) - rf^ | . (2) 

i—l ijj i—1 

Due to the isotopic character of the mixture, the interaction is the same between the 
different pairs of atoms. A simple representation of the potential is the Lennard- 
Jones interaction 



V{r) = 4e 



r / 



(- 

V r 



(3) 



where e = 10.22 K gives the depth of the potential and a = 2-556 A defines 
the length scale. Nowadays, the more accurate Aziz potential til and its revised 
version II2 HFD-B(HE) are used in realistic calculations. The He-He potential is 
characterized by a strong short range repulsion (such that the atoms, at a first 
approximation, can be considered as hard-spheres of diameter ~2.6 A) and a weak 
attraction at medium and large distances. An important feature to keep in mind 
is that the ionization energy of the atoms and the first excitation energy are large 
( of the order of the eV) when compared with the energies which play a role in the 
physics of Helium liquids and that are of the order of the Kelvin ( 1 eV ~ 11000 K). 
Therefore, the Helium atoms are truly the elementary constituents of the system. 
The chemical potential of the '^He impurity is defined by 

^ (^(3;A^4) I g(3;A^4) | ^(3;A^4)) _ (^(iV4) | \ ^{N^)) 

(vI/(3;A^4) I *(3;A^4)) (vI'(A^4) I *(iV4)) ' 

i.e. as the energy difference when one ^He atom is added to ''He atoms at 
constant volume. The quantities to be subtracted are of order A'4 while the result of 
order of unity. Therefore, one must explicitly take into account these cancellations 
between large quantities to have a good estimate of /ia . 

In a variational approach the next step consists in choosing a suitable trial 
function (^fp = ^'o(3;iV4)) for the ground state of A'4 "^He atoms plus one '^He in 
a volume 57, taken in the iV4, — > 00 limit, at constant ^He density, = N^/fl. 
The background wave function is simply obtained by omitting the impurity in the 
wave function. The extended Jastrow-Feenberg correlated wave function t2l, 

^'o(3;A^4) =i^2(3;A^4)i^3(3;iV4), (5) 
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represents a realistic choice for \l/o(3; -/V4). F2.3 are N-body correlation operators 
including explicit two- and three-body dynamical correlations. F2 is written as a 
product of two body Jastrow, ^He-^He and ''Hc-^'Hc correlation functions, 



(6) 



m>l=l,Ni 



and is given by a corresponding product of triplet correlations, 

The minimization of the ground state energy gives, in principle, the corre- 
lation functions via the solution of the Euler equations, (5£'o/(5/'"^^-' — and 
SEo/6f'^"''^''''^ = 0. Most of the results presented in this paper have been obtained 
by using an analytical form for the three-body correlation and solving the Eu- 
ler equations for the Jastrow factor within the Hypernetted Chain (HNC) integral 
equations technique for the distribution functions Eil. 

The parametrized triplet correlations have the Feynman form, 



/ 



(a,/3,7) 



{ra,rp,r^) = exp 



is 

eye 



^(r) being a parametrized variational function l3'E2I. An optimization of the triplet 
correlations in a larger functional space shows that the Feynman form is nearly 
optimal EJ. Optimal three body correlations for the impurity in the Ijiulli, solving 
the Euler equation for /("•'^■t), have been recently obtained and used EJu3. 

For the sake of illustration we give the energy equations with two body corre- 
lations only. The cancellations occuring in the calculation of may be effectively 
dealt with by writing the expectation value of i/(3; as Eq{3;N4) = i?! -|- /ig, 
where is the energy of the medium (proportional to N4) and /Z3 is the chem- 
ical potential of the ^He atom (of the order of unity). In fact, the background 
energy, E^, cancels in the difference (^. The background energy, obtained within 
the Jackson-Feenberg prescription for the kinetic energy, is given by 



El 



Pi 



V{r) 



2m4 



V2 ln/(4.4)(^) 



M3 = f^int + fJ-rea is dccomposcd into the interaction term fj,i- 



V{r) 



4m 



red 



-V2ln/(3.4)(^) 



(8) 



(9) 



"nired being the reduced mass nired = 'm3m4/{m3 + TO4), and firea, due to the 
medium rearrangement, 

t;2 



V{r) 



2m. 



.V2 ln/(4'4)(r) 



The distribution functions are defined by 
1 



pI 



Jdn\^{3:N4)\^ ' 



(10) 



(11) 
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P4 /df7|vI/(3;7V4)P ■ ^ ^ 

One can also introduce the corresponding static structure functions, via the 
Fourier transforms: 

= S^^'^^k) = 1 + Pi J dre'^ " [g'^^'^Hr) - l] . (13) 

and 

5(3) (fc) = 5(3^4) (fc) ^ P4J dre'^ '' [s'^^^^W - 1] . (14) 
In the fc — > Hmit 

Hm S'(3)(fc) = -(1 + a) (15) 

where a is the excess volume parameter, i.e. the relative increment in the molar 
volume when the impurity is added to the system, keeping the pressure constant. 
The mobility of the lighter mass ^He is larger and its presence tends to increase 
the pressure. As a consequence, the volume must be enlarged to keep the pressure 
constant. The experimental value of a at the ^He saturation density is 0.284. 

The problem of calculating the expectation values has so been translated into 
that of the evaluation of the distribution functions. Cluster expansion and integral 
equations provide a viable and effective tool to perform this task. To this aim, one 
introduces the function h = p — 1 and expands the distribution functions in powers 
of h. The terms of the expansions are diagrammatically classified and summed up 
to infinite orders by Hypernetted Chain integral equations, whose properties have 
been pedagpgically reviewed by Fabrocini and Fantoni, for both Bose and Fermi 
systems clcj. 

The impurity problem has been often analyzed within the Average Correlation 
Approximation (ACA) E3. ACA is obtained by taking the same dynamical corre- 
lation functions for the impurity and for the medium. In this case, (7(3>4) _ 
and 

^(^^^)(.r.) = ^4^, (16) 



therefore 

/^f.?^ = 2e(p4)+(^-l)i(^^(p4) (17) 

where e(p4) and t('')(p4) are the total and the kinetic '^He energies per particle, 
respectively, prea is then given in ACA by 

,,ACA _ P{Pi) N (^o\ 

Prea " 6(^4) (IS) 

Pipi) being the thermodynamical pressure. In this way the chemical potential is 
expressed as 

p^^^ = e{p,) + iW(p,) ^ ^,(,,) iW(p,). (19) 

P4 \m3 J \m3 J 
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The chemical potential of the impurity in ACA is the chemical potential of pure 
^He corrected by a kinetic energy term properly scaled to take into account the 
different mass of the impurity. Actually this phormula^an be easily generalized to 
wave functions containing also n-body correlations EjO. A recent difusion Monte 
Carlo calculationtl, which contains optimized n-body correlations, has given t^— 
14.32 K and fi^ — —7.27 K at saturation density. Employing these values, one gets 
fJ't'^^ = —2.58 K. If the optimization of the correlations affecting the impurity is 
performed, the chemical potential value is closer to the experimental one. 

An alternative trial wave function for Bose-systems has been recently proposed: 
the so called Shadow wave function methodcJ. It has been succesfully, used for 
microscopic calculations on liquid and |aolid ^He at zero temperature E3 and has 
been generalized to the impurity systemQ. In the SWF approach the wave function 
is given by 

^swFiR) = JdS Fswf{R,S) (20) 

where R = (ra, ri, tat) are the coordinates of the particles and S = (sa, si, ...sV) 
is a set of auxiliary variables, representing centers of motion of the particles and 
integrated over the whole space. The interparticle correlations are embedded in 

FswF = MR) X fpM rs - S3 I) X [] /p,(| f, - |) x >^,{S). (21) 

i=\ 

fp{R) = rit^i n - rs I) X Y\-^j fp{\ n - rj \) is a. Jastrow factor, as well as 

fs{S) = Ylfj^ifsil Si - S3 \) X ni<j /s(l Si ~ Sj I). "Sswr includes correlations 
at all orders (pair, triplet, ...) through the shadow correlations, after integrating 
over 5*. One of the main features of 'i'swF is the possibihty to describe both liquid 
and crystalline phases within the same functional form, respecting Bose symme- 
try and translational invariance. The expectation values between these trial wave 
functions are usually calculated by Monte Carlo techniques. The extension of ACA 
to the shadow variables approach, i.e. hjz. taking the correlation factors fp, fs and 
fps as recently optimized for pure ^He also for the impurity correlations, gives 
^3 = — 2.43 K at saturation density, close to the ACA value obtained with the 
DMC results as input. In the case of the SWF the three-, four- and n-body corre- 
lations are generated through the correlations with the shadow variables, while in 
the correlated basis theory they are explicitely introduced in the variational wave 
function and in the DMC method they are incorporated and/or optimized in the 
DMC algorithm for solving the many-body Schrodinger equation. 

The density dependence of the pure "^He energy and of the ^He impurity chem- 
ical potential is reasonably well described by CBF if the correlations are properly 
optimized and the cluster expansion diagrams are fully summed, with the inclusion 
of the elementary ones. CBF shows a good agreement with the results obtained 
within the two other approaches briefly discussed above. Even if the finite size 
of the simulation box introduces some uncertainties in the Monte Carlo meth- 
ods and some limitations in the study the long range behavior of the distribution 
functions, DMC provides in priiijdplc the most accurate results. In fact, by using 
the IIFD-B(HE) Aziz potential lla it reproduces the ^He experimental saturation 
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Figure 1. Pure ^He and impurity static structure functions at the *He saturation density. 



density po, as already mentioned it gives p!"^^''^\pQ) — —7.27(1) K and provides 
pf^^\pQ) = -2.79(25) K and a = 0.284(10). The kinetic and potential energies 
of the pure liquid "^He at po are t(po) = 14.32(5) K and v(po) = -21.59(5) K, 
respectively. 

Before leaving this subsection, we show in Figure 1 the static-structure functions 
S{k) and 5(3) (fc). They are the main inputs in the calculation of the excitation spec- 
trum and have been obtained at the experimental saturation density with the Aziz 
potential and solving the Euler-Lagrange equations for the two-body correlation 
functions, keeping the triplet correlations fixed Eil. S''^\k) has the proper behavior 
at fc ^ (^5|). The AC A cannot reproduce this limit because S^^)jj^{k) — S{k) — 1. 

2.1 The ^He impurity as a probe in liquid ^He 

The difference of expectations values in Eq.(^ shows that ^3, as given by the equa- 
tion, is not necessarily un upper bound to the true chemical potential. However, 
if we assume that the trial wave function is the exact wave function of liquid ^He, 
then p^^^ provides an exact upper bound to p*^^^ , 



Mr*(P4)<Mr*(P4)+(^-i)t(P4). 



(22) 



Since pl^^^ipi) and p'^^^*{pi) are experimentally known, the previous inequality 
can be used to establish a model independent lower bound to the kinetic energy of 
liquid ^He, 



t(p4) >t^^(p4) 



ma 

7714 — TTl3 



Mr*(P4)-Mr*(P4)]. (23) 
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At the same time, an upper bound to the potential energy per particle can be 
determined by 

v{p^) < v^^ipi) = e-f*(p4) - t^^'iPi) (24) 

where e'^^P*(/94) is the "^He energy per particle. At the experimental saturation 
density, /xf ^(po) = -7.17 K and ^ll^P{po) = -2.78 K, therefore t^^ipo) = 13.4 K 
and v^^{po) = —20.6 K. These bounds are basically respected by all the theoretical 
calculations performed with the Aziz potential. Also the extraction of the kinetic 
energy from the most recentinclastic neutron scattering data is consistent with the 
kinetic energy lower-bound EiJ. 



3 The excitation spectrum of the '^He impurity 

3.1 Variational calculation 

Let us focus now on the variational description of the impurity excited states. 
For pure ^He, a simple trial wave function for an excited ^atate characterized by a 
momentum q may be built in terms of Feynman phonons E3 

*q(iV4) =P4(q)*o(iV4), (25) 
where 'I'o(-^4) is the ''He ground state and 

P4(q)= e^''""' (26) 

is the "^He density fluctuation operator. The excitation energy, Lu{q), of ^'q(A'4) is 
(vl/q(7V4) \H~Eo I *q(A^4)) _ (*o(iV4) I P^P^iH - Eo) \ ^oiNi)) 



uiiq) = 

(■^r^(NA) I nt \H - Kn. n„] I -^niNi)) ^^^^ 



The first term of the right hand side vanishes when '^^{Ni) is the exact ground 
state or contains optimal two-body correlations. The second term leads to the well 
known Feynamn dispersion relation: 

c.(.)^c..(.) = ^, (28) 

where S(q) is the static structure function defined in Eg. dl^) . This relation gives the 
correct low momentum (g < OAA~^) linear behavior of the phonon-roton spectrum, 
ijj[q) = hqvs, where Vg is the speed of sound in liquid ''He {vs{po) ~ 238m/s), 
provided that 

lim S{q) ^ (29) 

In an analogous way, the impurity excitated stated is given by 

«'q(3;A^4) =P3(q)«'o(3;iV4), (30) 
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where p3(q) — e'^'^ '^^ is the impurity excitation operator. The excitation energy is 

^ {^^i3;Ni) I g(3;jV4) -go(3;iV4) | ^q(3;iV4)) _ V , . 

(vI/q(3;iV4) I *q(3;A^4)) 2m3' ^ ' 

and one obtains a dispersion relation corresponding to the free particle, with rrig = 

1713 . 

A better ansatz, which takes into account the backflow of ^He atoms around 
^He, is provided by 

*f^(3;7V4) =P3(q)i^B(q;3;^4)«'o(3;iV4) (32) 

with 

FB(q;3;iV4)= Yi exp [iq • (r, - r3)r/(r3,)]. (33) 

Backflow correlations do not change the binding energy of the impurity and only 
affect the excitation spectrum, which, however, remains parabolic: 

^BF{q) ^ T, [1 + ai + a2 + a3] , (34) 

where 



■013 

02 



f-^red 



ai=pij drg^^'^Hr) (^2r^{r) + ^W{r)^ , (35) 
P4 J drg^'^^Hr) (^r^irf + 1 [r^iv'irjf + 2Tj{r)rr,' {r}]^ , (36) 

pI J dr3idr325^^'^'^^ (^V3im2 + ^ [''31??31^32^32(r31 • ^32)^ + 2773i7732r32] 

(37) 

where g^^^^''^^ is a three body distribution function. A good choice for the variational 
function ri{r3i) has been found to be 



r/(r) = Ao exp 



(38) 



and the parameters ro, ujo and Aq are varied to find the minimum. With this 
variational ansatz, one gets a q-independent effective mass of mg/mj:.^ 1.7. A 
minimization through a full functional variation with respect to ?7(?'3i) c3 does not 
change this result. 

3.2 Correlated perturbative approach 

A systematic way to incorporate effects more complicated than the two-body back- 
flow is provided by the Correlated Basis Function perturbation theory. In such an 
approach, CBF is used to build richer excitations on top of ^'q by allowing the total 
momentum to be shared between the impurity and the phonons in the medium. The 
underlying idea consists in incorporating the non perturbative correlations directly 
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into the basis functions and then developing a perturbative expansion in this basis. 
UsuaUy, since the correlated states are a good approximation to the eigenstates of 
the hamiltonian, the first low orders of the perturbation theory are sufficient for 
getting reliable and realistic results. 

In the case of the impurity, the correlated basis is classified according to the 
number of Feynman phonons in each intermediate state, 

^-qiqL.qjS; A^4) = Psi^L - - . . - q„ )p4 (Qi ) • .p4 (q« ) *o (3 ; A^4). (39) 

The multiphonon states will mix in the perturbative process with the unperturbed 
state ^I'q = P3(q)^'o- The states ( |39| ) are not orthogonal among each other and ei- 
ther they must be properly normalized following the procedure described in Re£.(34) 
or the perturbative expansion must be carried on in a non orthogonal basis E3. 

The results presented in this work have been obtained including one phonon 
(OP) and two independent phonon (TIP) intermediate states and all the pertur- 
bative diagrams corresponding to successive rescatterings of the one phonon states 
(OPR). While the correlation factors are intended to care for the short range modi- 
fications of the ground state wave function due to the strongly repulsive interatomic 
potential, the basic physical effect induced by the perturbative corrections may be 
traced back to different types of backflow around both the impurity and the "^He 
atoms. 

In order to construct the CBF perturbative series, the unperturbed and the 
interaction Hamiltonians are defined via their matrix elements: 

_ (VI', I H I _ 

where Eg = Eq{3;N4) + £0(9) + <^^{q) is the eigenvalue of H for the state with 
momentum q. The diagonal matrix elements of Hj are zero by construction and, 
therefore, there are no first order perturbative corrections. 

The Brillouin-Wigner series for the perturbative correction AEg to E^, appro- 
priate for non-orthogonal states, is given by: 

AP {Hqj - EpNqj){Hjq - EpNjg) 

^^^-2^ e;^^ + ' 

where Eq is the final energy and the overlap matrix elements, Nqj = {'^q \ 'ij), 
take care of the non-orthogonalization. The series is first expanded around the 
correction to the medium energy, AEq = E4 — E^, and then resummed in such a 
way to cancel all the terms diverging in the A'4 — » 00 limit and originating from 
both the cluster expansion of the matrix elements and the perturbative expansion 
itself. This pmcedure has been devised for the ground-state eneifgy of an infinite 
Fermi system l3 and then generalized to the case of the impurity u. 

The difficulties related to the divergencies can be avoided by working in the 
orthogonalized scheme, in fact the n-phonon states may be Schmidt-orthogonalized 
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to states with a lower number of phonons. For instance, the orthogonahzed OP 
state reads as: 

^ ^ u75 ' ^^^> 

while the two-phonon state, 5'q;qiq2, may be orthogonalized in a similar way to 
^q;qi+q2 ^^'^ ^q:qi2- The orthogonalizaton makes the convergence of the series 
faster as the non orthogonalized states have large mutual overlaps 

The non diagonal matrix elements of the hamiltonian, H, can be easily eval- 
uated by assuming that the two- and three-body correlations are solutions of the 
corresponding Euler equations. Within this assumption, matrix elements involving 
n phonons may be expressed in terms of (n -I- l)-body structure functions. As an 
example, the matrix element between the OP state and the one without phonons, 
|q), is given by 

(q|i/|q;qi) - -[N^Siq^r'/^^^q- q.S^^^ {q,) , (44) 

where S{q) and S^^\q) are the two-body structure functions of cqs. ( [l3| ) and 
The diagonal matrix elements have the particularly simple form: 

(q;qi..q„|i/|q;qi..q„} = E'^ + eo{q) + ^ wpiqt) (45) 

1—1,71 

withi^o" = (*o|i?|*o)/(*o|*o)- 

Finally, the impurity excitation energy is e{q) — £0(9) + where, within 

the truncation we have used, 

Aeiq) ~ Aeop(g) + Aeripiq) + AeoPRiq) ■ (46) 

The different terms in (^) represent contributions from the corresponding inter- 
mediate states. 

The OP and TIP contributions to Ae(q) are shown in Figure 5 of Ref.(5). We 
stress that the CBF based perturbative approach requires the evaluation of two 
levels of diagrams, the perturbative and the cluster ones. We use Brillouin-Wigner 
perturbation theory, so the correction itself depends on e{q) and the series must be 
summed sclf-consistently. For instance, the OP contribution is solution of 

. n = V \i^\H-Eo-e{q)\q;qi)\' 
^^^'^^ 2^ e{q)-eoi\q-q,\)-WF{q,) 



qi 



(27r)3 V2m3y J ^ NiSiqi)e{q)-eo{\q~qi\)^WF{qi) 
and the OP effective mass, at g = 0, is 



'"3 



47r2p4 2m3 3 y ^ s{q)n^ 

2m3 2m4 



(48) 



The spectrum obtained by taking only OP states is very close to the LP one with 
an effective mass similar to that given by the variational calculation with backflow 
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correlations. Actually, in several papers it has been pointed out that second order 
perturbative expansion JKith OP states introduces two-body backflow correlations 
into the wave function ErEj. We find m'^CQP) = 1.8 m^, in good agreement with 
an analogous CBF treatment by SaarelacS (TO3 ~ 1.9 m^) and with the backflow 
variational calculationaj. 

The matrix elements involving TIP states (whose lenghty expressions are given 
in Ref.(5)), involve, as already anticipated, the two- and three-body structure func- 
tions, i.e. the Fourier transforms of the two- and three-body distribution functions, 
3^^H''i2) and (^('^^(ri, r2, ra). g^^^ may be evaluated within several approximations, 
the rpost common of which are the convolution (CA) and the superposition (KSA) 
ones 113. The CA correctly accounts for the sequential relation between g'^^^ and g^^^ 
and factorizes in momentum space, '5'pjj(qi, q2, qs) = S {qi) S {q2) S (q^) . The SA 

factorizes in r-space, 5^^^(ri, r2, rs) = g^^H^i2)ff^^H^i3)ff^^H^23), and adequately 
describes the short range region. The momentum space factorization property 
makes the CA more suited to our perturbative study. 

The two approximations give 'm'^{CA) = 1.6 TO3 and m'^lKSA) = 2.1 at 
q = 0, with 0£. and TIP states. A calculation including a four-body correction 
(the Abe term lB) in the three-body distributioii function and the OPR contribution 
provides mg = 2.2 7713 at saturation density □. Moreover, a key ingredient for a 
correct behavior of e{q) in the large q sector is a good description of the ''He roton. 
This requires the use of the superposition approximation that, on the other hand, 
is not correct in the ^He phonon region. The inclusion of backflow correlations in 
the correlated intermediate states gave an overall agreement with the experimental 
''He spectrum c3 but largely increased both the difficulty of the evaluation of the 
matrix elements and the uncertainty related to the use of high order distribution 
functions. In the impurity case a good compromise (termed CAl), which does not 
require a big computational effort, was found by using the CA and the experimental 
values of the ^He spectrum, ujexpt{q), in the energy denominators. 

In Figure 2, we show the impurity spectrum in CAl, together with the exper- 
imental ^He and ''He curves. The OPR terms are included and the LP and MLP 
fits to texptW) are shown. Since the branch of the dynamical response due to the 
excitations of the low concentration ■^He component in the Helium mixtures over- 
laps the collective ''He excitation at q > 1.7 A^'^ llJ, texpt{q) is not known in that 
region. 

The CBF spectrum is very close to texpt{q) up to its merging into the ^He 
dispersion relation. For the 7 parameter in the MLP parametrization, the theory 
gives ^{CBF) ^ 0.19 A^ and m'^{CBF) = 2.1 mg. If the spectrum is parametrized 
in terms of a momentum dependent effective mass, e{q) — q"^ / 2m^{q) , then we 
find in'^{q = 1.7 A^^) = 3.2 TO3, with an increase of ~ 50% respect to the q = 
value. 

Beyond q ^ 1.9 A^^, the impurity quasiparticle is no longer an excitation 
with a well defined energy and it can decay into *He excitations and acquire a 
finite lifetime, r. A finite r-value refiects a non zero imaginary part of the ^He 
complex on-shell self-energy, W{q) — 3I]((j, e(g)). W{q), amplified by a factor 4 
and computed with only OP intepnediate states, is shown in Figure 2. Our W{q) 
is close to the one found in Refs.HB. A numerical extrapolation of the computed 
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Figure 2. CBF/CAl (triangles), LP and MLP ^He single particle energies (in K). Stars and 
circles are the impurity and '*He experimental data, respectively. Black triangles are extrapolated 
CBF/CAl values. Black diamonds give the imaginary part of the impurity self-energy (in K). 
Momenta in 



^CBpiq) into the roton region does not show any evidences of a ^He roton-hke 

structure. 

Shadow wave functions have been used also for studying the excited states of 
the impurity. In this case, an excited state of momentum q is constructed by 
associating the momentum to the impurity shadow variables, 

*q(^) = j dSF{R, S)~5^. (49) 

Similarly, a ^He excited states is built as: 

^1{R) = j dSF{R,S)a^. (50) 

The momentum is carried on by the shadow variables with 

5q = e*'^■^^ aq = 5^e''>-^^ (51) 

where §^^3 are shadow variables modified by an explicit back-flow correlation 

Sj = Sj + ^(sj - si)X{sji). (52) 

The two shadow wave functions are eigenstates of the momentum operator, both 
with momentum q, and arc therefore orthogonal to the ground state. However, 
they are not orthogonal between them. They must be orthogonalized and a 2 x 2 
Hamiltonian must be diagonalized for each momentum in order to properly consider 
the collective and the single-particle impurity branches. This procedure has been 
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followed by Galli et aLQ by a Monte Carlo algorithm and the final spectrum is com- 
patible with a MLP spectrum with parameters = 2.06 and 7 = 0.0314 A^. 
The spectrum is close to the LP form and rnj = 1.74 7713 when there is no back- flow 
term in the shadow variables (A(rij) = 0). 

Other methods have confronted themselves with the impurity problem in the 
many-body arena. We first mention diffusion Monte Carlo, who has recently given 
an effective mass of m^/m^ = 2.2 □ at saturation density, at zero momentum. An- 
other one is based on a dynamic theory that allows for time dependent correlations, 
whose equations of motion are determined through a minimum action principle. 
Thisj-theory has been succesfully applied to study the excitation spectrum of pure 
■*He M. The method is very closely related to CBF and, in a simplifying hypotesis 
known as " uniform limit approximation" , one recovers the second order expression 
in the OP space. The resulting spectrum is well adjusted to the MLP form with 
mg = 2.09 ma and 7 = 0.114 at saturation density U. 



4 Conclusions 



We have shown in this contribution that a variational theory employing correlated 
wave functions is able to provide a good description of the ground state of one 
■^He impurity in liquid '^He. In addition, a perturbative expansion in a correlated 
basis may give a quantitative picture of the impurity excitation spectrum, pro- 
vided the basis considers correlated states with two independent phonons and one 
phonon rescattering diagrams, which play a non marginal role at large momenta. 
The agreement with other many-body methodologies, as shadow wave functions, 
diffusion Monte Carlo and time dependent correlations, is quite good. This gives 
confidence in the possibility of extending the CBF theory to other, less studied, as- 
pects of the physics of liquid Helium. To stay in the field of the mixtures, first CBF 



analyses of the ine 
mentum transfers 



neutron scattering cross sections, both at low and high mo- 
have revealed aLiL[£iicouraging semiquantitative agreement 



with the available experimental results Il3e3. This topic needs to be more carefully 
investigated, also within a close interaction with the experimental teams working, 
or willing to work, in this subject. 

Finally, to conclude this presentation, we rest our case. 
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